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SECTION  1 
INTRODUCTION 


In  recent  years  a  formulation  for  the  calculation  of  stochastic 
effects  on  electromagnetic  waves  due  to  structured  ionization  has  been 
developed  under  the  sponsorship  of  the  Defense  Nuclear  Agency  (DNA)  [Rino, 
1982;  Wittwer,  1979,  1982]  for  trans-ionospheric  links  at  radio  frequen¬ 
cies.  In  this  regime,  the  frequency  of  the  wave  is  always  substantially 
above  the  plasma  frequency  of  the  ionosphere  so  that  the  wave  propagates 
essentially  along  the  line  of  sight  path.  Approximations  which  make  use 
of  the  large  carrier  to  plasma  frequency  ratio  are  valid,  and  it  is  possi¬ 
ble  to  obtain  accurate  expressions  for  such  stochastic  quantities  as  the 
coherence  bandwidth  and  coherence  time  (or  their  Fourier  conjugates,  delay 
width  and  doppler  spread)  in  terms  of  a  fairly  small  number  of  physical 
parameters. 


At  high  frequencies  (HF)  the  situation  is  more  complicated.  In 
this  case  one  is  interested  in  signals  which  are  reflected  from  the  ioni¬ 
zation.  This  complicates  the  problem  in  that  the  path  followed  by  the 
wave  is  no  longer  a  straight  line,  but  is  instead  curved  as  the  wave  prog¬ 
resses  through  regions  of  changing  electron  density.  This  bending  of  the 
path  may  be  more  or  less  severe  depending  on,  for  instance,  the  launch 
elevation  angle  of  the  wave  or  the  steepness  of  the  electron  density  gra¬ 
dients  encountered  by  the  wave.  Since  the  wave  is  reflected,  there  may  be 
some  region  in  which  the  plasma  frequency  approaches  the  carrier  frequen¬ 
cy.  The  high  frequency  approximations  which  allowed  a  tractable  theory  in 
the  trans-ionospheric  case  may  thus  become  invalid  over  some  region  along 
the  path.  For  highly  oblique  paths  in  stratified  ionization,  however,  the 
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wave  carrier  frequency  remains  substantially  above  the  plasma  frequency 
even  at  the  reflection  point,  and  the  high  frequency  approximations  may 
still  be  valid. 

We  wish  to  explore  the  extent  to  which  the  trans-ionospheric 
theory  can  be  applied  to  such  internally  reflected  paths  at  HF.  This 
means  determining  whether  the  theory,  which  we  will  refer  to  as  the  DNA 
scintillation  model,  is  capable  of  predicting  the  appropriate  level  of 
scintillations  for  a  physically  reasonable  set  of  parameters.  Data  is  now 
available  to  make  such  a  test. 

Measurements  of  signal  strength  as  a  function  of  frequency  and 
delay  were  made  under  the  sponsorship  of  DNA  by  SRI  International  in  an 
experiment  called  the  HF  Channel  Probe.  The  measurements  we  will  consider 
were  made  over  a  one  hop  path  in  the  polar  region.  The  link,  located  in 
Greenland,  consisted  of  a  transmitter  at  Narssarssuaq  and  two  receivers  at 
Thule,  1918  km  to  the  north.  Details  on  the  experimental  setup  and  mea¬ 
surements  are  contained  in  Basler  et  al.,  [1985].  Their  measurements 
yield  information  on  the  degree  of  spreading  in  delay  and  doppler  frequen¬ 
cy  suffered  by  signals  at  HF  due  to  naturally  occurring  ionization  struc¬ 
ture  in  the  polar  region.  These  delay  spreads  and  doppler  spreads  can  be 
compared  to  the  predictions  of  the  DNA  scintillation  model  for  various 
values  of  the  model  parameters  as  a  test  of  the  applicability  of  the  for¬ 
mulation  at  HF.  To  confirm  the  model  would  require  independent  measure¬ 
ments  of  the  model  parameters  taken  simultaneously  with  the  measurements 
of  delay  spread  and  doppler  spread.  Since  this  data  does  not  exist,  we 
will  not  be  able  to  completely  validate  the  model,  but  enough  is  known 
about  the  model  parameters  that  any  gross  failure  of  the  model  should  be 
apparent . 

The  DNA  scintillation  model  is  based  on  integrations  taken  along 
the  mean  ray  path  of  quantities  involving  the  electron  concentration. 


1 

1 


Thus  one  must  know  the  electron  density  profile  along  the  path  followed  by 
the  ray.  In  conjunction  with  their  measurements,  the  HF  Channel  Probe 
experimentalists  made  oblique  soundings  over  the  link.  From  these  oblique 
ionograms  the  required  electron  density  profiles  may  be  inferred.  After 
first  reviewing  the  nature  of  oblique  ionograms  and  previous  methods  for 
determining  electron  density  profiles  from  them  (ionogram  scaling),  we 
present  in  Section  2  our  procedure  for  fitting  these  ionograms  using  a 
particular  functional  form  for  the  ionosphere,  the  quasi-parabolic  layer. 

Given  the  ionospheric  profile,  ray  trajectories  over  the  link 
can  be  numerically  constructed  and  the  scintillation  model  calculations 
performed.  We  have  incorporated  the  DNA  scintillation  model  calculations 
into  the  Jones-Stephenson  ray  tracing  code  for  this  purpose  [Jones  and 
Stephenson,  1975].  In  Section  3  we  will  briefly  discuss  the  relevant 
parts  of  the  DNA  scintillation  model  and  explore  the  model  dependence  on 
its  parameters.  From  this  we  will  gain  insight  into  the  relative  sensi¬ 
tivity  of  the  model  parameters,  and  this  information  will  be  used  to  guide 
us  to  an  appropriate  starting  point  for  testing  the  model  against  the  HF 
Channel  Probe  measurements.  In  particular,  we  will  find  that  the  model  is 
substantially  more  sensitive  to  one  of  the  parameters,  the  spectral  index, 
than  to  the  others  (although  two  other  model  parameters  will  be  found  to 
be  sufficiently  sensitive  to  be  important). 

In  Section  4  we  will  present  our  analysis  of  the  HF  Channel 
Probe  data  and  test  the  DNA  scintillation  model.  After  first  selecting 
several  data  sets  which  span  a  variety  of  propagation  conditions,  we  will 
extract  delay  and  doppler  spreads  from  the  measurements.  The  results  of 
the  parameter  sensitivity  investigation  performed  in  Section  3  will  moti¬ 
vate  an  initial  test  of  the  model  in  which  the  spectral  index  is  chosen  as 
the  primary  adjustable  parameter.  The  result  of  this  initial  test  will 
lead  to  speculation  about  which  of  the  remaining  parameters  should  be  made 
variable.  A  quite  successful  four  parameter  fit  to  the  data  will  finally 
be  performed. 
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The  results  of  our  analysis  will  be  summarized  in  Section  5.  We 
will  conclude  that  the  DNA  scintillation  model  is  capable  of  describing 
the  appropriate  level  of  stochastic  effects  observed  on  highly  oblique 
internal  paths  at  HF  for  reasonable  values  of  the  model  parameters.  This 
conclusion  is  based  on  our  result  that  good  fits  to  the  HF  Channel  Probe 
data  are  obtained  for  model  parameters  that  are  consistent  with  expected 
values  of  these  parameters,  and  in  some  cases  with  independent  (although 
not  simultaneous)  measurements  of  these  parameters. 


SECTION  2 
IONOGRAM  SCALING 


2.1  INTRODUCTION. 

The  DNA  scintillation  model  for  the  calculation  of  stochastic 
signal  structure  parameters  requires  that  integrations  be  taken  along  the 
mean  ray  path  of  quantities  which  involve  the  local  electron  density.  For 
this  reason  we  need  to  know  the  electron  density  profile  along  the  path 
followed  by  the  ray.  This  information  may  be  gleaned  to  some  degree  of 
approximation  from  the  oblique  ionograms  taken  by  the  HF  Channel  Probe. 

The  process  of  determining  electron  density  profiles  from  iono¬ 
spheric  soundings  is  known  as  ionogram  scaling.  In  this  section  we  will 
describe  vertical  and  oblique  ionograms  and  explain  their  appearance. 
Methods  of  scaling  oblique  ionograms  will  be  briefly  discussed  followed  by 
a  presentation  of  our  approach  to  this  problem.  Finally,  we  will  discuss 
a  method  for  handling  absolute  delay  uncertainties  associated  with  the  HF 
Channel  Probe  oblique  soundings. 

2.2  PROPERTIES  OF  VERTICAL  AND  OBLIQUE  IONOGRAMS. 

An  electromagnetic  wave  at  a  single  frequency  propagates  at  the 
phase  velocity  v^=  c/u  where  c  is  the  speed  of  light  in  vacuum  and  w  is 
the  index  of  refraction  of  the  medium.  Ignoring  collisions  and  any  exter¬ 
nal  magnetic  fields,  the  index  of  refraction  of  the  ionosphere  is  real  and 
is  given  by  w  =  (1  -  f^/f2)1/2.  Here  f  is  the  frequency  of  the  wave  and 
f  is  the  plasma  frequency,  f2  =  (c2r  /*)  N  ,  where  r  is  the  classical 


electron  radius  r  =  2.82  *  10" 15  m,  and  N  is  the  electron  number  densi- 
e  e 

ty.  A  wave  at  a  single  frequency  carries  no  information,  and  in  general 
signals  must  be  built  up  as  a  superposition  of  waves  at  a  number  of  fre¬ 
quencies.  While  each  frequency  component  of  the  signal  propagates  at  its 


phase  velocity,  the  signal  itself  will  propagate  at  the  group  velocity 
Vg=  dw/dk  (w=2*f,  k=2"A  where  A  is  the  wavelength  of  the  frequency  compo¬ 
nent).  We  will  be  primarily  concerned  with  signals  which  have  a  predomi¬ 


nant  frequency  component  known  as  the  carrier  frequency.  When  no  external 
magnetic  fields  are  present,  the  group  velocity  in  this  case  is  vg  =  cu, 
with  u  evaluated  at  the  carrier  frequency. 


An  ionospheric  sounding  is  a  measurement  of  the  time  of  travel 
from  transmitter  to  receiver  of  a  signal  at  a  given  carrier  frequency. 
This  transit  time,  or  delay,  is  given  by 

T  =  /£  -  1/  dS  [1  -V'/Z  ,  ID 
g  f2 

where  the  integral  is  over  the  path  followed  by  the  signal  (the  ray 
path).  An  ionogram  is  a  plot  of  these  delays  as  a  function  of  carrier 
frequency.  In  the  case  that  the  waves  are  propagated  vertically  into  the 
ionosphere  and  return  to  the  transmitter  location,  it  is  convenient  to 
define  the  virtual  height,  h'=  ct/2.  The  virtual  height  is  the  height 
implied  by  the  measured  delay  if  the  signal  had  propagated  at  the  speed  of 
light,  and  is  larger  than  the  true  height  h^  of  the  signal  reflection 
point. 

Consider  a  signal  propagated  vertically  from  the  earth  incident 
upon  a  single  ionospheric  layer  with  a  unique  maximum  in  electron  densi¬ 
ty.  As  a  signal  propagates  into  the  layer,  the  local  plasma  frequency 
grows  as  the  local  electron  density  gets  larger,  causing  the  index  of 
refraction  to  get  smaller.  The  group  velocity  is  reduced  and  the  delay  is 
correspondingly  increased.  The  wave  is  reflected  at  the  point  where  the 
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electron  density  is  such  that  the  plasma  frequency  equals  the  carrier  fre¬ 
quency  of  the  signal.  If  the  carrier  frequency  is  larger  than  the  plasma 
frequency  of  the  layer  maximum,  the  wave  will  penetrate  the  layer. 


A  synthesized  vertical  incidence  ionogram  is  shown  in  Figure 
1  (solid  curve)  as  a  plot  of  virtual  height  as  a  function  of  carrier 
frequency.  The  electron  density  profile  used  is  also  shown  (dotted  curve) 
plotted  as  a  function  of  plasma  frequency.  The  true  height  of  reflection 
at  a  given  carrier  frequency  is  the  height  of  the  equivalent  plasma  fre¬ 
quency  on  the  bottom  side  of  the  layer.  One  can  see  that  the  difference 
between  virtual  height  and  true  height  becomes  arbitrarily  large  near  the 
layer  maximum. 


In  the  case  of  a  monotonic  electron  density  profile,  Eq.  (1)  can 
be  inverted  to  yield  the  electron  density  profile,  given  the  virtual 
height  of  reflection  [Budden,  sec.  10.12], 


W  =  ®  ^  h'(f pSino)  da 


For  more  complicated  profiles  or  multiple  layers,  the  inversion  of  Eq. 
(1)  can  only  be  done  approximately.  In  these  cases  it  may  be  more 
convenient  to  assume  some  functional  form  for  the  electron  density  profile 
and  to  adjust  the  parameters  of  that  function  until  the  measured  virtual 
heiqhts  are  obtained.  Procedures  for  scaling  vertical  ionograms  are 
discussed  in  a  special  issue  of  Radio  Science  (Vol.  2,  No.  10,  October 
1967). 


We  now  turn  to  the  case  of  oblique  propagation,  in  particular 
to  the  bistatic  configuration  where  the  receiver  is  at  some  distance  D 
(ground  range)  from  the  transmitter.  Figure  2  shows  a  synthesized 
ionogram  corresponding  to  the  vertical  ionogram  of  Figure  1  for  the 
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Narssarssuaq  to  Thule  link  of  the  HF  Channel  Probe  measurements  (0  =  1918 
km).  The  group  path  is  defined  as  P'  =  ct  where  t  is  the  total  time  delay 
for  signal  propagation  from  transmitter  to  receiver.  The  horizontal  axis 
is  the  signal  carrier  frequency  and  has  been  labeled  oblique  frequency  to 
distinguish  it  from  from  the  carrier  frequency  of  the  equivalent  vertical 
ionogram. 

There  is  a  one-to-one  mapping  of  vertical  ionograms  to  oblique 
ionograms.  In  the  simplest  case  of  a  flat  earth  and  plane  stratified 
ionosphere,  three  theorems  suffice  to  define  this  mapping  [Davies,  sec. 
4.2].  The  first  is  a  direct  consequence  of  Snell's  Law  and  is  known  as 
the  Secant  Law, 

fOb  =  fV  Sec  *  •  <3> 

Here  $  is  the  complement  of  the  launch  elevation  angle.  This  law  states 
that  if  f  is  the  carrier  frequency  of  a  signal  reflected  at  some  true 
height  for  vertical  incidence,  then  a  higher  frequency  f  ^  will  be  reflec¬ 
ted  at  that  same  true  height  for  oblique  incidence.  The  second  is  Breit 
and  Tuve's  Theorem,  which  states  that  the  group  path  P*  for  ground-sky- 
ground  transmission  is  equal  to  the  length  of  the  triangular  path  defined 
by  the  launch  and  arrival  elevation  angles  (with  the  ground  range  D  as  the 
base  of  the  triangle).  The  third  is  Martyn's  Equivalent  Path  Theorem 
which  states  that  the  height  of  this  triangle  is  equal  to  the  virtual 
height  h'  of  a  vertically  incident  signal  at  the  equivalent  vertical  fre¬ 
quency  f  given  by  the  Secant  Law.  These  relationships  yield  the  mapping. 


These  expressions  must  be  modified  for  a  spherical  earth  and  ionosphere 
(which  is  the  case  for  Fig.  2),  but  Eq's  (4)  suffice  to  explain  the  gross 
features  of  oblique  ionograms. 

The  most  noticeable  feature  in  comparing  the  oblique  ionogram 
of  Figure  2  with  its  vertical  incidence  counterpart.  Figure  1,  is  that 
over  a  large  frequency  range  there  are  two  received  signals  for  a  given 
carrier  frequency.  This  is  a  general  property  of  single  layers  with  a  uni¬ 
que  maximum  electron  density.  The  signal  with  the  longer  delay  is  called 
the  high  ray,  the  other  the  low  ray.  These  rays  coalesce  at  a  frequency 
called  the  junction  frequency.  Figure  3  shows  the  ray  paths  for  the  high 
and  low  rays  at  10  MHz.  The  ionospheric  layer  is  also  shown,  with  the 
dashed  line  at  the  layer  maximum.  The  high  ray  is  seen  to  ride  the  crest 
of  the  layer.  A  slight  increase  in  launch  elevation  angle  would  cause 
that  ray  to  penetrate  the  ionosphere. 

2.3  METHODS  OF  SCALING  OBLIQUE  IONOGRAMS. 

Methods  for  scaling  oblique  ionograms  fall  basically  into  two 
categories,  those  which  involve  transforming  the  oblique  ionogram  to  its 
vertical  incidence  equivalent  and  scaling  that,  and  those  which  scale 
directly  from  the  oblique  trace.  Each  of  these  categories  can  be  further 
subdivided  as  to  whether  a  true  height  analysis  is  performed  (e.g.  Eq. 
(2))  or  whether  some  functional  form  is  used  in  a  fit  of  the  ionosphere. 


The  transformation  of  an  oblique  ionogram  to  its  vertical  inci¬ 
dence  equivalent  for  a  spherical  earth  and  spherically  stratified  iono¬ 
sphere  unfortunately  depends  not  only  on  geometrical  considerations,  but 
also  on  the  actual  refraction  experienced  by  the  ray,  hence  on  the  elec¬ 
tron  density  profile  itself.  One  may  ignore  these  refractive  effects  and 
carry  out  the  transformation  using  geometrical  considerations  only  by 
modifying  Eq.  (4b)  to  its  spherical  geometry  counterpart  (R  is  the 
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radius  of  the  earth), 

(£!)2  =  ?R(R+h ' )  [1  -  cos(JU]  +  h’2  .  (5) 

c.  cH 

Figure  4  compares  the  vertical  ionogram  of  Figure  1  with  the  result  of 
transforming  the  oblique  ionogram  of  Figure  2  using  Eq.  (5)  and  Eq.  (3) 
with  <(>  taken  as  the  apex  half  angle  of  the  triangular  path. 


There  is  a  significant  discrepancy  in  the  transformed  ionogram  as  one 
approaches  the  critical  frequency  of  the  layer.  Refractive  effects  can  be 
approximately  taken  into  account  by  modifying  Eq.  (4a)  to  [Smith,  (1939)], 
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fob  =  fv  K  sec  ♦  .  (7) 

where  the  correction  factor  k  is  largely  due  to  the  difference  in  angle 
of  incidence  between  different  ionospheric  heights  but  is  also  a  function 
of  the  electron  density  profile.  The  factor  <  typically  ranges  from  1  to 
about  1.2  .  Analytic  expressions  for  <  can  be  obtained  for  specific 
functional  forms  for  the  electron  density  profile  [Basler  and  Scott, 
(1973)].  The  two  curves  of  Figure  4  can  be  brought  into  substantial 
agreement  with  a  correction  factor  of  about  <=1.06  .  Thus  it  is  possible 
to  transform  an  oblique  ionogram  to  an  equivalent  vertical  incidence 
ionogram  to  a  fair  degree  of  accuracy.  Once  this  is  done  one  may  apply  a 
true  height  analysis  or  some  fitting  procedure  to  infer  an  associated 
electron  density  profile. 

As  we  remarked  above,  it  is  also  possible  to  carry  out  a  true 
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height  analysis  directly  from  the  oblique  trace  [Reilly,  (1985)], 
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Comparison  of  oblique  to  vertical  ionogram  conversion  using 
geometry  only  (dotted  curve)  with  the  true  vertical  incidence 
ionogram  (solid  curve). 


procedure  basically  involves  numerically  inverting  an  expression  for 
oblique  paths  analogous  to  Eq.  (1),  building  up  the  profile  from  the 
bottom.  To  perform  such  a  true  height  analysis,  whether  directly  from  the 
oblique  trace  or  from  the  equivalent  vertical  trace,  requires  that  the 
delays  are  known  fairly  precisely.  Ionograms  often  exhibit  substantial 
spreading  in  delay,  so  some  method  of  obtaining  an  appropriate  value  for 
delay  at  each  frequency  (i.e.,  fitting)  must  be  administered. 

We  have  found  that  the  most  convenient  method  for  scaling  the 
HF  Channel  Probe  soundings  is  to  parameterize  the  ionosphere  and  apply  a 
fitting  procedure  directly  to  the  oblique  ionogram  data.  Because  of  the 
large  spread  in  delay  typically  observed  in  these  measurements,  fitting 
must  necessarily  be  done  at  some  stage  in  the  analysis  anyway.  By 
choosing  a  convenient  functional  form  for  the  ionospheric  profile  and 
adjusting  the  parameters  to  obtain  a  best  fit  to  the  measured  oblique 
ionogram,  we  avoid  the  approximate  and  unnecessary  step  of  transforming  to 
equivalent  vertical  incidence. 

2.4  FITTING  THE  OBLIQUE  IONOGRAM  TO  A  QUAS I -PARABOL IC  LAYER 

The  ionosphere  typically  exhibits  two  layers  which  are  more  or 
less  distinct.  The  lower  layer,  or  E-layer,  has  its  maximum  electron  den¬ 
sity  at  a  height  of  about  110  km.  The  maximum  electron  density  of  the 
higher  F-layer  is  generally  much  larger  than  that  of  the  E-layer  and 
resides  in  the  range  of  200  to  400  km.  During  daylight  hours  the  F-layer 
sometimes  shows  a  smaller  enhancement  below  the  maximum  electron  density, 
known  as  the  F-layer.  In  this  case  the  main  maximum  is  denoted  the  F  - 

layer.  To  properly  fit  the  ionosphere  then,  two  or  three  parametric  lay¬ 
ers  should  be  used.  We  will  be  concerned  with  HF  Channel  Probe  scattering 
function  measurements  associated  with  reflection  from  the  F-layer.  These 
measurements  were  usually  made  at  carrier  frequencies  reasonably  close  to 
the  junction  frequency  for  the  F-layer,  and  in  this  region  the  signals  are 


only  minimally  affected  by  any  underlying  E-layer  ionization.  Furthermore, 
any  splitting  of  the  F-layer  is  sufficiently  indistinct  in  this  data  so  as 
to  warrant  only  a  single  layer  fit.  Thus  for  the  purpose  of  reproducing 
an  ionosphere  for  analyzing  the  scattering  function  measurements,  it  is 
adequate  to  fit  the  oblique  ionogram  data  to  a  single  layer  accurate  in 
the  nose  region  of  the  trace. 


Fitting  the  oblique  ionogram  data  requires  repetitive  calcula¬ 
tion  of  ray  paths  and  the  corresponding  delays  (or  group  paths).  For  this 
reason  it  is  convenient  to  choose  a  functional  form  for  the  layer  for 
which  analytic  expressions  for  group  path  and  ground  range  exist.  The 
quasi -parabol ic  layer  [Croft  and  Hoogasian,  (1968)]  provides  such  expres¬ 
sions  for  a  spherical  earth  and  spherically  stratified  ionosphere.  Written 
in  terms  of  plasma  frequency,  the  functional  form  of  the  layer  is 


[  J]2  (^)2} 


•  rh<  r  <  rm  ] 


elsewhere 


where  r  is  the  distance  from  the  center  of  the  earth  and  the  parameters 


fc  =  critical  frequency  of  the  layer 


rm  =  earth  centered  radius  of  the  maximum 


ym  =  semi -thickness  of  the  layer 


and  r.  =  r  -  y  .  We  introduce  the  substitute 


m  ■'m 


A  =  1  -  (!c)2  +  (IfV 
f  % 


(10a) 
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B  =  -2 r  (Ic^b)2  (10b) 

™  fy 
•'m 

C  =  (I§i)2  -  r02  cos2Bo  (10c) 

% 

cosy  =  (-^2.)  cosSo  (lOd) 

rb 

where  r0  is  the  radius  of  the  earth,  B0  is  the  launch  elevation  angle,  and 
Y  is  the  corresponding  angle  at  the  bottom  of  the  ionosphere.  The 
expression  for  ground  range  is  then 

D  =  2r„(Y-B0)  -  (_ _ lB2±flC] - }  (11) 

4C[$inY  +  &  +  — B — ]2 
rb  (2/C) 

and  the  group  path  is  given  by 

P'=  2rbsinY  -  2r0sine0  +  Z.  {-rbsinY  — in  - B2  -_4  AC - }.(12) 

^  4/A  [2Arb+B+2rb/AsinY]2 

For  a  given  carrier  frequency  and  set  of  layer  parameters  it  is  possible, 
using  Eq.  (11),  to  find  the  launch  elevation  angles  corresponding  to  the 
high  and  low  rays  which  give  a  desired  ground  range  D.  Once  these  angles 
are  known,  Eq.  (12)  yields  the  corresponding  delays.  Eq.  (11)  is  tran¬ 
scendental  in  0Q  ,  so  it  is  necessary  to  apply  a  searching  method  to  find 
these  launch  elevation  angles  for  fixed  ground  range.  It  is  therefore 
important  to  know  the  largest  launch  elevation  angle  Bmax  for  a  given 
carrier  frequency  beyond  which  rays  will  penetrate  the  layer.  This  angle 
can  be  obtained  by  solving  the  equation  for  ray  apogee  height. 
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Art  +  Brt  +  C  =  0  , 


for  80  using  the  expression  for  the  maximum  apogee  height. 


r  = 
t.max  2A 


One  approach  for  obtaining  a  synthesized  oblique  ionogram  from 
a  quasi-parabolic  layer  of  qiven  layer  parameters  is  to  first  find  the 
junction  frequency  fj  and  the  corresponding  launch  elevation  angle  8j  . 
High  and  low  ray  values  at  lower  frequencies  can  then  be  searched  for  at 


respectively  higher  and  lower  launch  elevations,  bounded  by  8max  above 


and  zero  below.  Figure  5  shows  a  plot  of  ground  range  as  a  function  of 
launch  elevation  angle  for  a  quasi-parabolic  layer,  plotted  parametrically 
with  carrier  frequency.  The  layer  parameters  in  this  case  are  f  =  4.5  MHz, 


h  =  300  km,  y  =  50  km  (h  is  the  height  of  the  layer  maximum  above  the 
m  m  m 


earth's  surface,  h  =  r  -  r0  ,  with  r0=  6359  km,  the  approximate  value  of 
mm 


the  earth's  radius  at  Greenland).  For  a  fixed  ground  range  D0  (e.g.  1918 
km  for  the  Narssarssuaq  to  Thule  link)  there  will  be  one  curve  whose  mini¬ 
mum  value  is  at  precisely  that  ground  range.  The  corresponding  frequency 
and  elevation  are  the  desired  fj  and  8j  .  Writing 


F(f,B)  =  D(f,B)  -  D0 


(15a) 


G(f ,3)  =  3D/36|f  , 


(15b) 


the  conditions  for  obtaining  fj  and  3j  are 


F(f  i,8  ■,)  =  0  =  G(f  ,,8,) 


Newton's  method  applied  to  two  functions  of  two  variables  provides  a 
convenient  way  of  numerically  solvinq  this  problem,  that  is,  solvinq 
iteratively  the  set  of  equations. 
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Fiqure  5.  Ground  ranqe  as  a  function  of  launch  elevation  plotted  para¬ 
metrically  with  carrier  frequency  for  a  quas i -parabol ic  layer 
with  parameters  fc  =  4.5  MHz,  hm  =  300  km,  ym  =  50  km.  The 
horizontal  dotted  line  is  at  D  =  1918  km  and  just  touches  the 
curve  fj  *  ll.?63  MHz  at  6j  =  14 .83  deq . 
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^new_  ^old  +  [FgG  -  F  Gg] 


[FfGe-  FpGf]  'old 


(17a) 


gnew_  3old  +  [F  Gf  -  FfG] 

[FfGp-  FpGf]  'old  ’ 


(17b) 


where  the  subscripts  refer  to  partial  differentiation. 


Having  found  the  junction  values,  a  search  can  be  made  at  lower 
frequencies  for  the  launch  angles  of  the  high  and  low  rays  by  stepping 
away  from  Bj,  cutting  the  step  size  down  each  time  the  calculated  ground 
range  exceeds  D0,  until  a  desired  accuracy  in  D0  is  obtained.  Eq.  (12) 
then  yields  the  group  path  for  these  frequencies  and  angles.  Figure  6 
shous  synthesized  oblique  ionograms  for  quasi -parabol  ic  layers,  showing 
the  variation  of  the  trace  as  each  of  the  layer  parameters  are  increased 


from  the  base  values  fc=4.5  MHz,  hm=300  km,  ym=50  km  (solid  curve). 


The  HF  Channel  Probe  oblique  ionogram  dala  taken  by  SRI  Interna¬ 
tional  consist  of  power  measurements  at  128  delays  for  each  of  576  carrier 
frequencies.  One  such  measurement  is  displayed  in  Figure  7.  This  sound¬ 
ing  was  made  on  March  20,  1985  (day  79  refers  to  the  day  number  of  year 
1985)  at  18:31  GMT,  corresponding  to  a  local  time  of  approximately  16:00. 
Only  data  within  about  10  dB  of  the  peak  measured  power  is  shown.  A  fair¬ 
ly  well  defined  F-layer  is  apparent  as  well  as  some  evidence  for  an  under¬ 
lying  E-layer.  The  data  in  the  upper  left  corner  of  the  plot  is  the  two 
hop  echo  of  the  F-layer.  The  F-layer  trace  exhibits  a  slight  nose  exten¬ 
sion.  This  is  a  scatter  effect  frequently  observed  in  high  latitudes. 


We  fit  the  measured  ionogram  data  to  a  quasi-parabol ic  layer 
utilizing  a  nonlinear  regression  fitting  routine.  This  routine  synthe¬ 
sizes  an  oblique  ionogram  correspond  ng  to  a  set  of  starting  layer  param- 
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Figure  6.  Oblique  ionogram  dependence  on  quasi -parabol ic  layer  parameters 
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eters,  and  a  power  weighted  chi  square  to  the  delay  data  is  calculated 
(the  weighting  is  done  in  dB).  The  layer  parameters  are  then  changed 
slightly  and  the  direction  toward  decreasing  chi  square  in  parameter  space 
is  determined.  After  stepping  the  parameters  in  this  direction,  the  pro¬ 
cess  is  repeated  until  chi  square  no  longer  decreases. 

The  result  of  one  such  nonlinear  regression  fit  to  the  data  of 
Figure  7  is  shown  in  Figure  8  .  Signals  corresponding  to  the  two  hop  mode 
were  ignored  in  the  fit  by  using  only  data  between  9.5  and  13  MHz. 
Similarly,  the  E-layer  data  was  windowed  out  by  using  only  data  with  a 
group  path  greater  than  2050  km.  The  data  in  this  fit  has  power  greater 
than  -30  dB  with  a  peak  power  near  -20  dB  (relative  to  some  standard).  The 
fit  minimizes 

/  p.  measured  _  p,  QP  . 

X2  =  l  Weight.  ( - 1—— - -1 - -)  (18) 

1  Error 

where  Weighti  =  Power.  -  Threshold,  Threshold  =  -30  dB,  and  the  error 
associated  with  each  delay  was  arbitrarily  set  to  10  km  (0.033  msec). 
Since  all  data  points  are  qiven  the  same  error,  any  other  constant  error 
would  yield  the  same  fit.  Contributions  to  x2  for  data  beyond  the 
junction  frequency  of  the  quasi-parabolic  layer  were  calculated  using  the 
group  path  of  the  junction  frequency  for  P'w  .  Thus  the  fit  is  not  overly 
affected  by  nose  extensions. 

2.5  HANDLING  ABSOLUTE  DELAY  UNCERTAINTIES. 

We  will  be  examining  data  from  the  first  two  HF  Channel  Probe 
campaigns,  in  particular  data  taken  on  October  3,  1984  (Day  277,  fall 
campaign)  and  March  20,  1985  (Day  79,  spring  campaign).  All  data  from 
these  two  campaigns  are  subject  to  an  approximate  200  usee  uncertainty  in 
absolute  time  delay  measurements  due  to  timing  system  drifts  in  the 


315  MHz 


experimental  apparatus  as  well  as  uncertainties  associated  with  the 
Transit  satellite  system  used  to  update  the  system's  internal  clocks 
[Basler,  Roy  P.,  private  communication].  This  corresponds  to  a  60  km 
uncertainty  in  group  path.  Quasi-parabolic  fits  to  the  ionosphere  can 
vary  substantially  over  this  range.  As  an  example  we  show  in  Table  1  the 
results  for  quasi-parabol  ic  fits  of  the  ionogram  shown  in  Figure  9a  over 
the  range  At  -  ±  200  usee. 

Table  1.  Quasi-parabolic  fits  for  forced  absolute  delay  shifts  At  of  the 
measured  data.  The  data  and  At  =  0  fit  are  shown  in  Fig.  9a.  The  last 
column,  column,  x2/N.  .  is  a  measure  of  the  goodness  of  the  fit  (see  Eq. 
(18)).  Ndata  is  the  number  of  data  points  used  in  the  fit. 


At  (usee) 

fc  (MHz) 

h  (km) 
m 

Xn  (kn,) 

*2/Ndata 

-  200 

3.98 

207.6 

28.7 

11.95 

0 

4.47 

278.5 

22.5 

9.82 

+  200 

5.03 

356.8 

21.8 

20.35 

Such  differences  produce  substantial  changes  in  ray  path,  and 
as  we  shall  see  in  Section  3,  quantities  calculated  in  the  DNA  scintilla¬ 
tion  model  can  scale  quite  strongly  with  path  length  and  electron  densi¬ 
ty.  An  estimate  of  absolute  delay  can  often  be  made  by  using  the  E-layer 
trace  by  assuming  that  the  bottom  of  the  trace  lies  at  a  height  of,  say, 
100  km.  However  the  E-layer  trace  is  not  always  visible  in  the  ionogram 
data.  We  note  from  Table  1,  though,  that  the  fit  obtained  in  the  At  =  0 
case  was  significantly  better  than  for  the  other  two  cases,  as  is  evi¬ 
denced  by  the  x2  of  the  fit.  Thus  there  is  reason  to  believe  that  the 
delay  shift  is  relatively  small  in  this  example. 

It  is  possible  to  include  absolute  delay  shift  as  a  parameter  in 
the  fitting  process.  Figure  9b  shows  the  fit  obtained  for  the  preceding 


25 


*»  V! 


example  when  this  was  done.  Note  that  positive  Ax  corresponds  to  addi¬ 
tional  delay  added  to  the  delay  data  and  is  not  equivalent  to  applying  an 
opposite  shift  to  the  layer  height  parameter.  The  delay  shift  of  at  = 
13.3  usee  (4  km  in  group  path)  is  indeed  small,  but  enough  to  improve  the 


fit  by  five  percent  to  a  chi  square  of  X  'N<jata  =  9.35. 


Table  2  summarizes  results  from  this  kind  of  analysis  for  two 
sequences  of  ionograms.  The  corresponding  fits  are  shown  in  Figures  9 
through  14.  In  all  cases  the  fits  were  improved  by  including  the  delay 
shift  (as  is  evidenced  by  comparison  of  the  resultant  chi  squares), 
although  the  improvement  is  small  in  some  cases.  An  independent  check  on 
the  success  of  this  fitting  procedure  in  obtaining  appropriate  delay 
shifts  is  to  check  the  resulting  position  of  the  E-layer  trace.  The  first 
sequence  (Day  277)  was  chosen  as  an  example  where  no  E-layer  trace  appears 
in  the  thresholded  data  (recall  that  a  threshold  of  -30  dB  was  used  for 
fitting  purposes).  Indeed,  examination  of  the  unthresholded  ionograms 
(not  shown)  shows  only  faint  E-layer  traces.  These  traces  are  nearly  cor¬ 
rectly  positioned  in  the  17:31  and  17:46  GMT  cases  (assuming  that  the 
bottom  of  the  E-layer  is  at  100  km),  but  the  18:01  GMT  E-layer  trace  is 
slightly  higher  than  the  previous  two  cases.  The  delay  shifts  obtained  in 
the  Ax  fit  of  this  sequence  are  sufficient  to  produce  consistency  in  these 
E-layer  traces,  even  though  no  E-layer  information  was  used  in  the  fitting 
procedure.  In  the  second  sequence  (Day  79),  some  evidence  for  the  E-layer 
appears  near  the  bottom  of  the  plots.  A  virtual  height  of  100  km  corres¬ 
ponds  via  Eq.  (5)  to  a  group  path  of  1942  km,  and  while  the  bottom  of  the 
E-layer  trace  is  not  really  discernible  in  these  plots  because  of  thresh¬ 
olding,  comparison  with  the  At  =  0  cases  shows  that  a  large  improvement  in 
the  positioning  of  the  E-layer  has  been  made  by  fitting  the  delay  shift 
(even  in  the  Day  79  18:01  case,  which  showed  the  least  improvement  in  x2). 
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Table  2.  Comparison  of  ionogram  fits  for  At  =  0  and  At  varied. 


Day/time 

fc(MHz) 

h  (km) 
m 

ym(k») 

At  (usee) 

X2/Nd 

277  17:31 

4.47 

278.5 

22.5 

0 

9.82 

17:46 

4.72 

306.7 

22.1 

0 

9.61 

18:01 

4.79 

329.4 

21.8 

0 

10.90 

17:31 

4.48 

280.4 

22.1 

13.3 

9.35 

17:46 

4.71 

305.3 

21.3 

16.4 

9.36 

18:01 

4.74 

313.7 

21.9 

-44.4 

10.12 

079  18:01 

5.29 

360.2 

45.3 

0 

8.92 

18:16 

5.34 

345.3 

33.8 

0 

9.25 

18:31 

5.32 

348.0 

38.2 

0 

8.11 

18:01 

5.10 

333.4 

49.1 

-80.2 

8.83 

18:16 

5.10 

311.7 

37.0 

-87.1 

7.18 

18:31 

5.08 

317.7 

37.6 

-88.6 

7.12 
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Figure  11.  Comparison  of  ionogram  fits  for  no  delay  shift  Ar 
and  Ar  fitted  (b) . 
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OBLIQUE  IONOGRAM  Day  79  18:31:  0  (hr,min,sec) 


fc  =  5.315  MHz 
hm  =  348.0  km 
ym  =  38.2  km 
A  r  =  0.0  msec 
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ym  =  37.6  km 
A  t  =  -0.0886  msec 
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Figure  14.  Comparison  of  ionogram  fits  for  no  delay  shift  Al¬ 
and  Ar  fitted  (b) . 
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Of  the  six  cases  examined  in  the  two  sequences  above,  three  will 
be  used  in  our  analysis  presented  in  Section  4  (Day  277  17:46,  Day  79 
18:16,  Day  79  18:31).  Three  other  cases  not  shown  here  will  also  be 
analyzed  in  Section  4  (see  Figures  26,  28,  and  32).  The  appearance  of 
E-layer  traces  in  the  first  two  of  these  cases  leaves  little  doubt  that 
even  though  large  delay  shifts  are  obtained,  they  are  appropriate.  The 
last  case  shows  only  a  hint  of  an  E-layer  trace  in  the  full  unthresholded 
ionogram  (not  shown),  but  the  resultant  delay  shift  seems  to  be  consistent 
with  positioning  the  bottom  of  this  faint  trace  at  1942  km  group  path. 
However,  we  note  that  delay  and  doppler  spread  data  associated  with  this 
trace  is  most  poorly  fit  in  the  analysis  of  Section  4,  and  thus  the  layer 
parameters  may  be  suspect.  Even  so,  we  conclude  overall  that  the  best  fit 
ionospheric  layer  obtained  by  incorporating  At  as  a  fit  parameter  is  prob¬ 
ably  the  appropriately  positioned  ionospheric  layer. 
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SECTION  3 

STOCHASTIC  SIGNAL  STRUCTURE  PARAMETERS 


3.1  INTRODUCTION. 

The  high  latitude  ionosphere  is  characterized  by  structured 
ionization  aligned  along  the  earth's  magnetic  field.  These  structures,  or 
striations,  are  thought  to  be  the  result  of  structured  electron  precipita¬ 
tions  which  produce  large  scale  irregularities  that  are  subsequently  con- 
vected  through  the  auroral  zone  and  polar  caps.  Smaller  scale  structures 
are  produced  as  plasma  instabilities  (e.g.  E*B  or  gradient  drift  instabil¬ 
ities)  break  up  the  large  scale  structure  during  transport  [Kelley  et  al., 
1982],  Rapid  expansion  along  the  field  lines  accounts  for  the  quasi- 
cylindrical  shape  of  the  striations.  Structure  scale  sizes  parallel  to 
the  field  may  be  more  than  ten  times  as  large  as  the  perpendicular  scale 
sizes.  It  is  the  structure  perpendicular  to  the  field  which  dominates  in 
producing  scintillations  of  radio  signals.  These  amplitude  and  phase 
fluctuations  are  the  result  of  distortions  imposed  on  the  wave  front  as  it 
propagates  through  these  irregularities. 

The  DNA  scintillation  model  was  formulated  as  a  relatively  sim¬ 
ple  parametric  description  of  the  statistical  effects  of  scatter  from  such 
structure.  In  this  section  we  review  the  stochastic  signal  structure 
parameters  relevant  to  the  HF  Channel  Probe  data  and  their  calculation  as 
formulated  by  Wittwer  [Wittwer,  1979,1982].  The  dependence  of  the  model 
on  its  parameters  will  be  mapped  out,  providing  information  on  the  sensi¬ 
tivity  of  the  stochastic  parameters  to  the  model  parameters.  This  will 
give  insight  into  which  p.irameters  may  be  safely  fixed  and  which  should  be 
varied  when  we  fit  the  HF  Channel  Probe  data  in  Section  4. 
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3.2 


THE  SCATTERING  FUNCTION. 


The  distortions  imposed  on  the  wave  front  by  structured  ioniza¬ 
tion  may  be  viewed  in  the  ray  optics  approximation  as  micromultipath. 
Energy  reaching  the  receiver  in  a  particular  mode  (e.g.  F-layer  high  ray 
or  low  ray)  includes  contributions  from  many  paths,  each  the  result  of  a 
sliqhtly  different  sequence  of  scatterings  from  the  structure.  Signal 
components  traveling  over  different  paths  will  in  general  have  different 
propagation  times,  and  thus  the  received  signal  will  be  spread  in  time  of 
arrival  or  delay.  Temporal  variation  of  the  channel  structure  results  in 
a  corresponding  frequency  (doppler)  broadening  of  the  signal.  Because  the 
irregularity  structure  of  the  ionosphere  must  be  described  statistically, 
a  statistical  treatment  of  the  delay  and  doppler  spreading  is  required. 
The  DNA  scintillation  model  provides  such  a  description.  The  HF  Channel 
probe  has  measured  the  received  power  spectrum  in  delay  and  doppler  for 
transmitted  pseudo  random  noise.  This  spectrum,  known  as  the  scattering 
function,  provides  a  direct  measure  of  the  spreading  in  delay  and  doppler 
due  to  the  channel,  as  we  will  see  below. 

The  response  of  a  channel  to  a  transmitted  signal  Eo(t)  may  be 
modeled  linearly  [Bello,  (1963)], 

ao 

E(t)  =  /  dx  h( t ; t )E0 ( t-x)  .  (19) 


The  function  h(x;t)  describes  the  response  of  the  channel  at  time  t  to 
a  signal  transmitted  at  the  earlier  time  t-x.  If  the  transmitted  signal 
were  an  impulse,  E0(t)  =  <5  ( t ) ,  then  h(x;t)  would  itself  be  the  received 
field.  Hence  it  is  called  the  time-variant  impulse  response  function.  It 
depends  on  time  because  of  the  temporal  variation  of  the  channel.  The 
dependence  on  delay  is  a  reflection  of  multipath.  The  Fourier  transform 
of  the  impulse  response  function,  the  time-variant  channel  transfer 
function,  describes  the  carrier  frequency  response  of  the  channel. 


2TTTT 


(20) 


H(f;t)  =  /  dT  h( t ;t )  e_i 


It  represents  the  response  of  the  channel  to  a  single  frequency 
sinusoid. 

To  obtain  a  statistical  description  of  the  channel  response  we 
consider  correlations  of  these  channel  system  functions.  Consider  the  two 
frequency  -  two  time  channel  correlation  function, 

r(f 1 »f2  ;ti ,t2 )  =  <H* ( f 1  ;t ! )  H(f2;t2)> 

=  /  dT  7  d€  <h*(T;ti)  h(C;t2)>  el2ir(f  lT"f ^  .  (21) 


Two  useful  approximations  can  be  made.  The  first  is  that  signals  arriving 
with  different  delays  (different  paths)  are  uncorrelated.  This  uncorrel¬ 
ated  scatter  (US)  approximation  takes  the  form 

<h*( t ; t ! )  h( C;t2)>  =  *h( T;t ! ,t2)  <s( t-c)  .  (22) 

As  a  result,  the  two-frequency  -  two-time  correlation  function  becomes  a 
function  only  of  the  frequency  separation  Af  =  f  x  -  f2. 


r(^i»ir2JtiJt2)  "  r ( ; 1 1 , 1 2 ) 


=  /  dt  *h(T;t1,t2)e 


i  2  TT  Af  T 


(23) 


The  second  approximation  is  to  assume  that  the  statistics  of  the  channel 
do  not  change  with  time  so  that  the  correlation  functions  depend  only  on 
time  separation. 
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*h(T;ti,t2)  =  *h(T,At) 


(24) 


where  At  =  ti  -  t2.  If  this  approximation  is  valid  the  channel  is  said  to 
be  wide  sense  stationary  (WSS).  Wide  sense  stationarity  allows  the  time 
dependence  of  the  correlation  functions  to  be  Fourier  analyzed  in  a  single 
doppler  frequency  domain  f  d .  We  define  the  scattering  function  S(x,fd)  by 

00  i  2  tt  f  ,  At 

*h(i,At)  =  /  dfrf  S(T,fd)  e  ,  (25) 

.00 


and  thereby  obtain  under  the  assumption  of  wide  sense  stationary  - 
uncorrelated  scatter  (WSSUS)  the  double  Fourier  transform  relationship. 


7  7  /  i2"Afx  i2nf  .At 

r ( Af ; At )  =  J  dx  J  dfd  S( t , f d )  e  e  d 


Since  r(Af;At)  is  a  correlation  function  in  both  of  its  argu¬ 
ments,  S(x,fd)  is  properly  interpreted  as  a  power  spectrum  in  both  argu¬ 
ments.  To  see  this  we  consider  the  autocorrelation  of  the  received  field 
at  two  different  times. 


OO  00 

<E*(t)  E(t-At)>  =  /  dx  /  d5  <h*(x;t)  hU;t-At)>  <E0*(t-x)E0(t-At-5)> 

—  00  _  OO 


00  OO 

=  /  dx  /  dc  4>h(x;At)  <5(x-5)  <E0*(t-x)E0(t- At-C)>  (27) 

-OO  -00 


=  /j dx  4>h(x,At)  <E0*(t-x)E0(t-At-x)> 


where  in  the  first  line  we  have  assumed  that  the  statistics  of  the  chan¬ 


nel  are  independent  of  the  transmitted  signal  and  in  the  second  line  that 
the  WSSUS  assumption  holds.  In  terms  of  the  scattering  function,  this  is 


★  "  “  i2uf  .At  * 

<E  (t)  E(t-At)>  =  /  dx  /  df  ,  S ( t , f , )  e  <E0  ( t- t)  E0(t- At-x>.( 28) 

_oo  .»  oo 


Setting  At  =  0  yields  the  expression  for  the  received  signal  strenqth  in 
terms  of  the  transmitted  signal  strength. 


<| E( t ) | 2>  =  /  dx  /  dfrf  S(x,fd)  <| E0( t- t) | 2> 


The  scattering  function  is  thus  the  del ay-doppler  power  spectrum  of  the 
received  signal  due  to  the  channel. 


In  practice,  the  scattering  function  may  be  measured  by  trans¬ 
mitting  a  wideband  signal  and  cross  correlating  the  received  signal  with  a 
delayed  version  of  the  transmitted  signal.  For  transmitted  random  noise 
this  yields  the  impulse  response  function  (up  to  an  amplitude  factor  that 
can  be  included  in  the  definition  of  h(x;t)). 


<E*(t)  E0 ( t - t ) >  =  /  dC  h*(S;t)  <E0*(t-5)E0(t-T)> 


=  /  dS  h  U;t)  6(x-0 


h  (  x ;t ) 


Correlating  this  result  at  separated  times  yields  t  (x,At)  if  WSSUS  is 


assumed.  Finally,  Fourier  transforming  4>.  over  At  to  the  doppler  domain 


gives  the  scattering  function  S(x,f  ). 
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A  sample  scattering  function  measurement  from  the  HF  Channel 
Probe  is  shown  in  Figure  15.  The  peak  power  in  this  plot  has  been  nor¬ 
malized  to  unity.  The  delay  origin  is  at  6.7  msec  and  has  been  adjusted 
to  allow  for  a  200  usee  delay  introduced  by  the  experimental  apparatus. 
Two  distinct  modes  are  apparent.  Comparison  with  the  ionogram  of  Figure 
10,  taken  one  minute  earlier,  shows  that  these  modes  correspond  to  the  low 

and  high  rays  of  the  F-layer.  Each  mode  exhibits  a  substantial  spread  in 

both  delay  and  doppler,  characterized  by  standard  deviations  which  we 
denote  respectively  and  o^.  The  doppler  spread,  is  conventionally 
defined  as  the  1/e  point  of  the  scattering  function  (i.e.,  for  a  gaussian 
=  /2c<d).  Because  of  the  Fourier  transform  relationship  between  the 
scattering  function  and  the  two  freguency  -  two  time  channel  correlation 
function.  Eg.  (26),  these  widths  imply  correspondi ng  coherence  widths  in 

the  conjugate  variables,  i.e.  for  gaussian  scattering  functions, 

f0  =  l/(2nox)  tq  =  l/(irvd)  .  (31) 

The  coherence  bandwidth  f0  represents  the  frequency  spread  over  which  sig¬ 
nals  will  remain  correlated.  Similarly,  the  decorrelation  time  t0  is  a 
measure  of  the  time  over  which  signals  will  remain  correlated.  Both  of 
these  quantities  are  of  fundamental  importance  in  the  design  and  perform¬ 
ance  assessment  of  HF  communication  systems.  Broad  delay  and  doppler 
spreads  translate  to  narrow  coherence  bandwidths  and  short  decorrelation 
times.  This  situation  is  likely  to  occur  under  nuclear  conditions,  and  in 
order  to  design  systems  which  will  be  free  of  intersymbol  interference  in 
such  cases  one  must  be  able  to  estimate  these  signal  structure  parameters 
given  information  about  the  environment.  The  DNA  scintillation  model  was 
formulated  to  provide  such  a  character i zat ion  for  transionospheric  satel¬ 
lite  communication  links.  The  HF  Channel  Probe  scattering  function 
measurements  afford  us  the  opportunity  to  test  the  appl icabi 1 ity  of  this 
formulation  on  internally  reflected  paths  at  HF  over  a  variety  of  condi¬ 
tions. 
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3.3 


THE  DNA  SCINTILLATION  MODEL. 


In  order  to  obtain  a  statistical  description  of  signal  propaga¬ 
tion  effects,  one  must  have  a  statistical  description  of  the  environment 
affecting  the  propagation.  It  is  therefore  necessary  to  consider  the 
correlation  function  of  electron  density  variations  in  the  ionosphere, 

r AN ( ^ )  =  <ANe^  ANe^  +  ^  (32) 

where  ANg(r)  =  Ng(r)  -  <Ng(r)>.  The  three-dimensional  Fourier  transform 
of  this  correlation  function, 

*N(k)  =  / d3P  elk  *  p  rAN(p)  (33) 

is  a  spectral  density  function  that  is  readily  related  to  the 
three-dimensional  power  spectral  density  (PSD)  of  index  of  refraction 
fluctuations.  The  connection  is  made  by  considering  the  approximation 

u  =  [1  -  fp2/f2]1/2  -  1  -  ( 1/2 ) ( f p2/ f 2 )  (34) 

with  f  2  =  (c2r  /ti )  N  (r  is  the  classical  electron  radius,  r  =  2.82  * 
p  e  e  e  p 

10"15m).  A  deviation  in  electron  density  ANg  produces  a  corresponding 
deviation  in  the  index  of  refraction, 

AU  =  -  C(A2re)/(2*)]  ANe  .  (35) 

Thus  the  PSD  for  index  of  refraction  fluctuations  is  related  to  that  of 
the  electron  density  variations  by. 


PSD(k)  =  [(X2re)/(2ir)]2  *N(ic) 


(36) 


The  DNA  scintillation  model  is  based  on  a  power  law  representation  of  this 
PSD,  which  for  axi symmetric  striations  may  be  written  as 

8*3/2<am2>l^2  L|(  r(n)/r(n-3/2) 

PSD(i<a,k1)  =  (1  +  1_^2  k^2  +  L(|2  k#2jn  .  (37) 

In  this  expression,  is  the  spatial  wave  vector  perpendicular  to  the 
earth's  magnetic  field,  k^  is  the  spatial  wave  number  parallel  to  the  mag¬ 
netic  field,  Lx  and  L|(  are  the  structure  outer  scale  sizes  relative  to  the 
magnetic  field,  <Am2>  is  the  index  of  refraction  variance,  and  n  is  the 
power  law  index.  The  corresponding  one-dimensional,  or  in  situ  PSD  var- 
ies  as  l/(x  ^)  for  large  k,  so  the  spectral  index  (SI)  is  defined  to  be 
SI  =  2n-2  .  The  normalization  of  Eq.  (37)  requires  n  >  1.5  ,  and 
Wittwer's  formulation  [Wittwer,  1979,  1982]  is  valid  in  the  range  1.5  <  n 

£  4- 

The  value  of  the  spectral  index  depends  both  on  the  striation 
profile  and  on  the  distribution  of  striation  sizes  [Wortman  and  Kilb] 
While  a  value  of  SI  =  2  has  been  advocated  for  a  number  of  years  on  the 
basis  of  several  in  situ  measurements  (see  references  cited  in  Wittwer, 
1979),  more  recent  DNA  HILAT  measurements  have  yielded  values  (averaged 
over  more  than  a  year  in  the  region  of  the  HF  Channel  Probe  experiment)  of 
SI  =  1.5  (day),  1.63  (night)  with  a  significant  probability  of  occurrence 
throughout  the  range  0.5  to  2.5  [Basler  et  al.,  1985]. 


While  Wittwer's  formulation  of  the  DNA  scintillation  model 
allows  the  calculation  of  a  number  of  stochastic  signal  structure  param¬ 
eters  (Rayleigh  phase  variance,  mean  square  log  amplitude  fluctuation,  S4 
scintillation  index,  parallel  decorrelation  time,  angle  of  arrival  vari- 


ance),  we  will  concentrate  on  the  two  relevant  to  the  HF  Channel  Probe 
scattering  function  measurements  the  delay  spread  o  and  the  decorrelation 
time  perpendicular  to  the  ray  path,  which  we  label  x0  (actually  t0  =  min 
but  for  the  paths  we  consider  the  perpendicular  decorrelation 
time  is  always  smaller  than  the  parallel  decorrelation  time).  Again,  the 
doppler  spread  is  related  to  the  decorrel ation  time  by  vd  =  l/(irr0). 


A  differential  phase  variance  can  be  obtained  from  Eq.  (27)  by 
considering  correlations  of  the  phase  shift  imparted  in  traversing  seg¬ 
ments  of  the  scattering  region  larger  than  the  correlation  distance  of  the 
index  of  refraction  fluctuations  [Wittwer,  1979].  This  then  assumes  that 
the  ray  remains  in  the  scattering  region  for  distances  longer  than  this 
correlation  length  in  the  ray  direction.  Over  the  Narssarssuaq  to  Thule 
path,  the  magnetic  dip  is  near  80  deg.  and  the  angle  between  the  ray  path 
and  the  magnetic  field  lines  typically  will  not  be  less  than  60  deg.  in 
the  refracting  region.  Thus  we  are  far  from  being  parallel  to  the  magnet¬ 
ic  field  lines  where  such  an  assumption  might  not  be  valid.  Letting  ds 
represent  the  differential  geometric  ray  pathl&ngth,  the  differential 
phase  variance  is 


8?5/2T(n-l)  re  ^<N  >2  l-l  t  (38) 

r(n-3/2)  k2  e  (cos2'!' +  a02sin2ip) x/2 


where  ao  =  L (| /L^  and  is  the  angle  between  the  ray  direction  and  the  mag¬ 
netic  field  lines.  We  have  assumed  that  the  index  of  refraction  variance 
is  related  to  the  electron  density  variance  by  <Au2>  =  (2*re/k2)2  <ANe2>. 
This  approximation  assumes  that  the  carrier  frequency  of  the  signal  is 
always  sufficiently  larger  than  the  local  plasma  frequency  so  that  the 
approximation  y  =  (l  -  f^/f2)1/2  »  1  -  ( 1/? ) ( f p  2/f 2 )  is  valid.  This  will 
be  true  even  near  the  reflection  point  for  the  highly  oblique  paths  we  are 


considering.  We  have  also  assumed  that  the  electron  density  variance  is 
proportional  to  the  square  of  the  average  electron  density. 


will  not  repeat  it  here,  but  refer  the  reader  to  the  references  [Wittwer, 
1979,  1982].  The  constant  Bn  is  given  by 


IM-  cl  ♦  hoi 

if1/2r(n-l)  3  (n-!. 


where 


h(e)  = 


ln(e) 


n  =  2 


— 1 — , C 1  -  el4-2n|]  otherwise  , 
4-2n|  J 


(n-l)“° • 14  n  <  2 


( n  - 1 ) 


-0  30 


n  >  2 


c  -  e'1'1  Kn-l(c) 
2n'2  r(n-l) 


and  e  is  the  ratio  of  inner  to  outer  structure  scale  sizes.  For  e  small, 
the  hyperbolic  Bessel  function  may  be  expanded. 


Kp_i(e)  -  (1/2)  r(n-l)  (c/2) 


-( n- 1 ) 


which  yields  cp  ■  1  . 


The  signal  decorrelation  time  is  estimated  in  the  frozen  plasma 
approximation  (assuming  fixed  transmitter  and  receiver)  as 


B(n,a|)  Lj_  ^  j  t  ^  da$  +  cos2^  +  ao2sin2^jm'2 ^ -  1/m 


cos2')'  +  ao2sin2ij> 


where  vx  is  the  velocity  of  plasma  motion  perpendicular  to  the  ray  path, 
and  which  we  have  taken  to  be  constant  and  removed  from  the  integral.  The 
function  B(n,a^)  helps  insure  the  appropriate  behavior  for  small  phase 
variances,  and  is  given  by 


B(n,o*)  =  min[l,  (-0.34  n2  +  2.51  n  -  2.00)  ( cr^B  )1/m] 
<p  q>  n 


The  parameters  required  for  the  calculation  of  delay  and  doppler 
spreads  in  this  formulation  are  then  n  (or  equivalently  the  spectral  index 
SI  =  2n-2) ,  l_x,  L n ,  vx,  £  =  <ANg2>/<Ne>2,  and  e.  Of  these,  vx  affects 


ly  the  calculation  of  doppler  spread  (or  t0),  and  enters  this  linearly. 


PARAMETER  DEPENDENCE. 


As  we  have  just  seen,  six  parameters  are  required  by  the  DNA 
scintillation  model  to  calculate  delay  and  doppler  spreads.  Testing  the 
model  against  the  HF  Channel  Probe  data  requires  adjusting  these  param¬ 
eters  in  a  fit  of  the  measured  scattering  functions.  A  single  mode  yields 
only  two  data  points,  a  delay  spread  and  a  doppler  spread,  and  hence  only 
two  parameters  can  be  fit  unless  some  relationship  between  the  parameters 
of  different  modes  is  assumed.  This  possibility  will  be  discussed  further 
in  Section  4,  but  for  the  present  we  will  assume  that  each  mode  will  be 
fit  independently.  Since  vx  appears  only  in  the  doppler  calculation,  dop¬ 
pler  spreads  can  always  be  fit  by  adjusting  vx  after  the  delay  spread  has 
been  fit.  Thus  of  the  remaining  five  parameters,  four  must  be  fixed  and 
one  varied  in  fitting  the  delay  data.  For  this  reason  it  is  important  to 
discern  which  of  these  five  parameters  is  the  most  sensitive. 


We  have  mapped  out  the  parameter  dependence  of  the  DNA  scintil¬ 
lation  model  for  a  nominal  ionosphere,  that  corresponding  to  the  oblique 
ionogram  of  Figure  2.  This  quasi-parabolic  layer  has  a  critical  frequency 
of  4.5  MHz  at  300  km,  and  a  half  width  of  50  km.  The  ray  paths  for  the  10 
MHz  rays  used  were  shown  in  Figure  3.  The  high  ray  travels  a  significant 
distance  at  nearly  the  peak  of  the  layer,  and  scintillations  calculated 
over  the  high  ray  will  be  correspondingly  larger  for  the  same  parameters 
than  for  the  low  ray.  Figures  16  -  25  show  the  variation  of  delay  spread 
and  doppler  spread  with  each  of  the  parameters.  Base  values  for  the 
parameters  were  SI  =  2,  =  10  km,  L,  =  100  km,  5  =  0.2,  e  =  0.001,  and 

v^  =  50  m/sec. 

Figures  16  and  17  show  a  very  strong  dependence  on  the  spectral 
index  in  the  vicinity  of  SI  =  2.  For  example,  the  delay  spread  varies  by 
an  order  of  magnitude  for  the  low  ray  and  two  orders  of  magnitude  for  the 
high  ray  for  SI  ranging  from  1.5  to  2.5  .  The  dependence  on  ,  shown 
in  Figures  18  and  19,  is  also  fairly  strong.  However  for  /L(|  <*  0.05  - 

0.1,  an  expected  reasonable  range,  the  variation's  less  than  a  factor  of 
two  for  any  of  the  quantities.  Figures  20  and  21  show  almost  no  depend¬ 
ence  on  Ln  once  it  is  two  or  three  times  larger  than  L^.  Thus  L(|  may  be 
safely  fixed  at  100  km  without  affecting  the  fit  of  the  HF  Channel  Probe 
data.  The  ratio  of  the  inner  to  outer  structure  scale  size,  e,  is  poorly 
known,  but  Figures  92  and  23  show  the  dependence  on  this  parameter  to  be 
exceedingly  small  over  many  orders  of  magnitude.  It  enters  into  the  cal¬ 
culation  of  the  coefficient  Bn  logarithmically  for  n=2  (the  case  shown 
here)  and  as  a  low  power  otherwise,  but  in  any  case  the  dependence  is 
weak.  We  will  fix  e  at  0.001. 

The  linear  dependence  of  the  delay  spread  and  the  square  root 
dependence  of  the  doppler  spread  on  S,  shown  in  Figures  24  and  25  is 
easily  understood.  From  Eq.  (38)  we  note  that  do^/ds  scales  as  5  <Ng>2. 
Examination  of  Eq.  (10)  then  shows  that  the  doppler  spread  will  scale  as 
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Figure  17.  Delay  and  doppler  spread  variation  with  spectral  index  for 
the  high  ray  of  Figure  3.  The  remaining  parameters  are 
fixed  at  Lj_  =  10  km,  l_n  =  100  km,  E  =  0.2,  e  =  0.001, 

Vi  =  50  m/sec. 
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ad  -  /s  /£  <Ng>  for  n  =  2.  The  scaling  of  delay  spread  can  be  similarly 
obtained  by  examining  Eq.  (40).  While  it  was  not  explicitly  shown,  the 
Rayleigh  phase  variance  scales  as  ~  s  £  <Ng>2  (i.e.  it  goes  like  the 
integrated  phase  variance).  The  second  term  in  Eq.  (40)  scales  as  s4  ?2 
<Ng>4  for  n  =  2.  Remembering  to  take  the  square  root  of  Eq.  (40)  to  ob¬ 
tain  the  delay  spread.  Figures  24  and  25  clearly  show  that  it  is  this 
second  term  which  dominates.  While  these  scaling  rules  change  slightly  as 
n  becomes  less  than  2,  it  is  clear  that  the  dependence  on  £  is  sufficient¬ 
ly  slowly  varying  to  be  of  secondary  concern  in  parameter  fitting. 

These  figures  show  that  the  spectral  index  is  substantially  the 
most  sensitive  parameter  in  the  formulation.  Given  the  uncertainty  in  our 
present  understanding  of  an  appropriate  value  for  the  spectral  index,  or 
more  correctly,  the  range  of  values  over  which  this  parameter  has  been 
observed  to  vary,  it  makes  most  sense  at  this  stage  to  vary  the  spectral 
index  in  fitting  the  HF  Channel  Probe  data  while  fixing  the  other  param¬ 
eters  to  reasonable  values.  or  £  might  be  varied  secondarily  when  com¬ 
paring  data  from  associated  modes. 

The  scaling  rules  just  obtained,  especially  the  one  for  delay 
spread,  ~  s2  K  <Ne>2,  show  that  the  formulation  is  very  sensitive  to 
the  ray  path,  in  particular  for  the  high  ray  which  may  have  a  significant 
path  length  near  the  layer  critical  frequency.  As  we  saw  in  Section  2.5, 
the  uncertainty  in  the  ionogram  delay  measurements  can  result  in  large 
uncertainties  in  the  fitted  layer  parameters.  Hence  there  is  some  danger 
in  fitting  the  scattering  function  measurements  that  ray  path  inaccuracies 
may  be  absorbed  into  the  fitted  parameters. 


SECTION  4 

TESTS  OF  THE  SCINTILLATION  MODEL 


4.1  INTRODUCTION. 

In  the  previous  sections  we  have  presented  a  method  for  scaling 
the  HF  Channel  Probe  oblique  ionograms  and  have  described  the  DNA  scintil¬ 
lation  model  calculation  for  delay  and  doppler  spreads.  We  now  test  this 
model  against  the  HF  Channel  Probe  scattering  function  measurements.  Our 
approach  will  be  to  trace  rays  through  the  quasi -parabol ic  ionosphere  ob¬ 
tained  from  an  ionogram  fit,  and,  with  the  scintillation  model,  calculate 
delay  and  doppler  widths  over  these  paths.  We  seek  a  set  of  model  param¬ 
eters  which  reproduce  the  measured  scattering  function  widths. 

This  section  begins  with  a  description  of  how  delay  and  doppler 
spreads  are  extracted  from  the  data.  We  then  consider  initial  fits  to  the 
low  ray  data  and  obtain  some  insight  into  how  to  handle  associated  low  and 
hiqh  ray  modes.  Finally,  we  consider  fits  of  these  associated  modes  and 
discuss  the  significance  of  the  resulting  model  parameter  values. 

4.2  THE  DATA. 

The  fitting  methods  used  in  our  analysis  force  us  at  present  to 
look  at  each  case  individually.  This  makes  the  analysis  sufficiently 
tedious  that  it  is  impossible  to  examine  the  whole  volume  of  data  supplied 
by  the  HF  Channel  Probe  experiment.  Hence  we  have  selected  six  data  sets 
which  represent  a  variety  of  propagation  conditions  to  analyze  in  detail. 


An  example  scattering  function  measurement  was  shown  in  Figure 
15.  This  measurement  clearly  displayed  two  distinct  modes,  associated 
with  the  low  and  high  rays  reflecting  from  the  F-layer.  Splitting  between 
ordinary  and  extraordinary  modes  is  not  evident  in  this  example  and  is 
only  suggested  in  a  few  of  the  cases  we  have  examined,  so  we  confine 
ourselves  to  the  examination  of  ordinary  waves  only  (i.e.,  we  ignore  the 
external  magnetic  field  of  the  earth  when  calculating  ray  trajectories, 
although  it  will  be  included  in  the  calculation  of  stochastic  effects). 

The  shape  of  the  scattering  function  suggests  that  it  may  be 
well  represented  by  a  sum  of  two  bivariate  gaussians.  Assuming  such  a 
functional  form  for  the  scattering  function  provides  a  convenient  way  to 
obtain  delay  and  doppler  widths  from  the  measurements.  We  assume  the 
form. 


where  the  subscripts  1  and  2  refer  to  the  low  and  high  ray  modes  respec¬ 
tively.  We  fit  (least  squares)  the  measured  scattering  function  data  to 

this  form  to  obtain  the  best  fit  parameters  P.,  u..,  v  u  o  ..  Of 

i  di  di  ti  ti 

these,  we  are  primarily  interested  in  the  widths  v^.  and  a^.  which  we  take 
to  be  the  measured  doppler  and  delay  spreads. 

The  selected  data  sets,  gaussian  fits,  and  corresponding  iono- 
grams  are  shown  in  Figures  26  through  37  and  Tables  3  through  8.  Note 
that  different  mesh  sizes  are  used  in  the  fitted  scattering  function  plots 
compared  to  the  measured  scattering  function  plots  (half  the  frequency 
resolution  and  twice  the  delay  resolution).  The  powers  are  normalized  to 
a  peak  power  of  one.  Comparison  of  the  scattering  functions  with  the 
associated  ionograms  shows  in  each  case  that  the  modes  are  properly 
interpreted  as  low  and  high  ray  F-layer  reflections. 
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Figure  27.  Fit  to  two  bivariate  gaussians  of  the  scattering  function 

measurement  of  Figure  26.  Final  fit  parameters  are  shown  in 
Table  3. 

Table  3.  Bivariate  gaussian  fit  parameters  for  the  scattering  function 
of  Figure  26,  Day  277,  13:47  GMT. 
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Figure  29.  Fit  to  two  bivariate  gaussians  of  the  scattering  function 

measurement  of  Figure  28.  Final  fit  parameters  are  shown  in 
Table  4. 

Table  4.  Bivariate  gaussian  fit  parameters  for  the  scattering  function 
of  Figure  28,  Day  277,  15:18  GMT. 
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Figure  31.  Fit  to  two  bivariate  gaussians  of  the  scattering  function 

measurement  of  Figure  30.  Final  fit  parameters  are  shown  ir. 
Table  5. 

Table  5.  Bivariate  gaussian  fit  parameters  for  the  scattering  functior 
of  Figure  30,  Day  277,  17:47  GMT. 
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Mqure  35.  pit  to  two  bivariate  qaussians  of  the  scattering  function 

measurement  of  Fiqure  34.  Final  fit  parameters  are  shown  in 
Table  7. 

Table  7.  Bivariate  qaussian  fit  parameters  for  the  scattering  function 
of  Figure  34,  Ray  79,  1 B : 19  GMT. 
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Figure  37.  Fit  to  two  bivariate  gaussians  of  the  scattering  function 
measurement  of  Figure  36.  Final  fit  parameters  are  shown 
in  Table  8. 

Table  8.  Bivariate  gaussian  fit 
of  Figure  36,  Oay  79, 

parameters  for  the  scattering  function 
18:32  GMT. 

P;  -  1- 

P2  -  0.9351 

H  -  0.5113  Hz 

di 

i/J  -  -0.6392  Hz 

di 

-  1.778  Hz 

di 

-  1.073  Hz 

d2 

u  -  6.924  msec 

ri 

v  ,  -  7 . 393  msec 

d2 

a  -  33.94  jisec 

a  -  21.62  /jsec 

4.3 


INITIAL  FITS. 


The  DNA  scintillation  model  parameter  dependence  investigation 
performed  in  Section  3.4  suggested  that  the  spectral  index  is  the  most 
sensitive  parameter  in  the  model.  As  a  starting  point  then,  we  have  fit 
exactly  the  low  ray  delay  spread  by  fixing  the  other  parameters  to  reason¬ 
able  values  and  adjusting  the  spectral  index.  The  doppler  spread  for  the 
low  ray  was  then  exactly  fit  by  adjusting  the  plasma  velocity  parameter. 
These  resulting  parameters  were  subsequently  used  to  calculate  delay  and 
doppler  spreads  for  the  high  ray.  Comparison  of  the  calculated  high  ray 
results  to  the  data  will  be  used  to  gain  insight  into  a  better  fitting 
method  for  these  associated  modes  to  be  performed  in  the  next  subsection. 

We  have  fixed  the  structure  outer  scale  sizes  to  10  km,  L  = 
100  km,  the  ratio  of  inner  to  outer  structure  sale  size  to  e  =  0.001,  and 
the  fractional  electron  density  variance  to  £  =  0.2  .  Table  9  lists  the 
spectral  index  SI  and  the  plasma  velocity  vd  required  to  fit  each  of  the 
low  ray  modes  exactly  as  well  as  the  results  obtained  when  these  param¬ 
eters  are  applied  to  the  high  ray  (listed  in  the  table  as  the  ratio 
[calculated  spread]/data) .  The  resulting  spectral  index  values  are 
somewhat  large  compared  to  the  ONA  HILAT  results  quoted  in  Section  3.3  (SI 
w  1.5  to  1.63),  but  are  still  well  within  the  observed  bounds  of  the  HILAT 
measurements.  The  plasma  velocities  are  also  very  reasonable  when 
compared  to  electron  drift  velocity  measurements  made  usinq  the  Sondre 
Stromfjord  incoherent-scatter  radar,  located  on  the  path  at  about  one 
third  of  the  link  distance  from  the  HF  Channel  Probe  transmitter  (Basler 
et  al.,  1985].  As  yet  we  have  been  unable  to  obtain  the  velocity  data 
from  the  days  of  these  scattering  function  measurements  (although  it 
exists),  but  we  can  say  that  velocities  ranqinq  from  eero  to  400  m/s  (and 
often  much  higher)  are  typically  observed. 
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Table  9.  Result  of  adjusting  the  spectral  index  and  plasma  velocity 
parameters  to  fit  the  low  ray  modes  exactly  and  subsequent 
application  of  these  parameters  to  the  high  ray  modes 
(listed  as  the  ratio  [calculated  spread]/data) .  Fixed 
parameters  are  Lx=  10  km,  LB=  100  km,  e  =  0.001,  5  =  0.2. 
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Somewhat  disturbing,  however,  is  the  observation  from  Table  9 
that  delay  and  doppler  spreads  calculated  for  the  high  ray  using  these 
parameters  (which  fit  the  low  ray  quantities  exactly)  are  in  every  case 
considerably  too  large.  These  results  are  shown  graphically  in  Figure 
38.  Here  we  plot  the  calculated  doppler  spread  divided  by  the  measured 
doppler  spread  as  a  function  of  the  calculated  delay  spread  divided  by  the 
measured  delay  spread  (circles).  The  solid  curve  is  the  square  root  of 
the  delay  axis.  We  see  that  in  five  of  the  six  cases  the  doppler  spread 
discrepancy  is  quite  close  to  the  square  root  of  the  delay  spread  discrep¬ 
ancy.  We  note  from  the  analysis  of  Section  3.4  and  the  scaling  rules 
obtained  there  that  there  is  one  parameter  that  scales  in  precisely  this 
way,  the  fractional  electron  density  variance,  S.  Figure  38  suggests  that 
allowing  this  one  parameter  to  vary  between  the  low  and  high  ray  trajec¬ 
tories  will,  to  a  substantial  degree,  bring  the  calculated  values  of  both 
delay  and  doppler  spread  into  agreement  with  the  data. 

4.4  FITTING  ASSOCIATED  MOOES. 

The  parameters  of  the  power  spectral  density  of  index  of  refrac¬ 
tion  variations  (PSD,  Eg.  37)  are  considered  locally  h  mogeneous  quanti¬ 
ties  [Tatarskii,  1961,1971],  which  means  that  they  may  be  slowly  varying 
functions  in  any  spatial  direction.  In  fact,  there  is  no  a  priori  reason 
to  expect  any  of  the  model  parameters  to  be  essentially  constant  over  any 
particular  ray  path,  much  less  over  different  ray  paths.  The  results  of 
the  initial  tits  described  in  the  previous  subsection  suggest  that  to  qet 
agreement  between  model  calculations  and  measured  data  at  least  one  of  the 
parameters  will  have  to  be  considered  variable.  We  have  suggested  that 
the  fractional  electron  density  variance  f,  may  be  an  appropriate  choice, 
althouqh  one  miqht  do  As  well  by  varyinq  the  spectral  index  or  L(  (since 
one  can  argue  that  the  ft w  points  in  Fig  ire  38  may  be  fit  with  any  number 
of  curves).  Tince  our  model  ionosphere  is  spherically  stratified,  the 


vertical  direction  is  the  only  direction  in  which  such  variations  can  be 
justified  in  the  model.  However,  the  peak  altitudes  for  the  low  and  high 
ray  trajectories  differ  by  only  10  to  25  km  in  the  cases  we  are  examin¬ 
ing.  These  distances  are  a  fraction  of  the  parallel  outer  scale  size, 
which  is  the  appropriate  scale  size  to  consider  for  vertical  variations 
since  the  magnetic  di  is  so  high  in  the  polar  region.  The  spectral  index 
should  be  only  slowly  varying  over  distances  defined  by  the  outer  scale 
size  in  that  direction  for  the  formalism  to  be  valid,  and  the  same  is  true 
about  the  scale  sizes  themselves.  The  expression  assumed  for  electron 
density  variances  (Eq.  39),  on  the  other  hand,  is  a  low  order  approxi¬ 
mation  designed  to  give  a  rough  order  of  magnitude  account  of  the  layer  as 
a  whole.  It  should  not  be  expected  to  hold  in  a  detailed  fashion  over 
different  modes.  A  more  accurate  model  for  these  variances  may  be  expect¬ 
ed  to  include  some  spatial  dependence  in  £  (perhaps  scaling  as  the  gradi¬ 
ent  of  the  mean  electron  density). 

As  a  first  cut  at  introducing  a  spatially  dependent  £,  we  con¬ 
sider  separate,  but  constant,  values  for  the  high  and  low  rays.  This 
should  be  roughly  equivalent  to  using  a  smoothly  varying  £  since  each  ray 
path  within  the  layer  can  be  well  represented  by  a  constant  altitude  seg¬ 
ment  (see  e.g..  Figure  3).  All  other  parameters  will  be  considered 
constant  over  the  entire  layer.  Simultaneously  fitting  associated  low  and 
hiqh  ray  modes  allows  four  parameters  to  be  varied  since  there  are  four 
data  points  (two  delay  spreads  and  two  doppler  spreads).  We  choose  to  fit 
the  spectral  index,  plasma  velocity,  and  the  fractional  electron  density 
variances  £x  and  £2  (low  and  high  rays  respectively). 

The  results  of  these  fits  are  shown  in  Table  10.  Fractional 
electron  density  variances  for  the  high  ray  modes  overall  are  an  order  of 
magnitude  smaller  than  those  for  the  low  ray.  This,  though,  corresponds 
to  only  a  factor  of  three  or  four  in  the  root  mean  square  quantity 
<>Ne/Ne.  The  plasma  velocities  are  little  changed  from  the  values  obtained 


Table  10.  Results  of  a  four  parameter  simultaneous  fit  of  the  low  and 
high  ray  modes.  Fit  parameters  are  the  spectral  index  SI, 


plasma  velocity 
variance  for  the 
parameters  are  L 
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low  and  high  rays,  £1  and  £2- 
=  10  km,  L „=  100  km,  e  =  0.001. 
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in  our  initial  fits  of  the  previous  subsection.  The  spectral  index  val¬ 
ues,  however,  have  changed  significantly  and  are  now  very  close  to  the 
mean  values  inferred  from  the  DNA  HILAT  measurements.  It  is  gratifying  to 
see  that  the  freedom  introduced  by  making  S  a  fit  parameter  has  worked  in 
such  a  way  as  to  move  the  spectral  index  in  this  direction. 


SECTION  5 
CONCLUSIONS 


We  have  tested  the  ONA  scintillation  model  calculations  of  delay 
width  and  doppler  spread  against  the  HF  Channel  Probe  measurements  of 
these  quantities  for  the  case  of  high  frequency  waves  propagated  through 
the  naturally  occurring  ionization  structure  of  the  polar  region.  Our 
analysis  has  required  the  scaling  of  oblique  ionograms  taken  over  the 
channel  probe  link  in  order  that  a  model  ionosphere  could  be  constructed 
through  which  rays  could  be  numerically  traced  and  the  model  calculations 
performed.  Measured  delay  widths  and  doppler  spreads  were  extracted  from 
the  data  by  fitting  the  measured  scattering  functions  to  two  bivariate 
gaussians  corresponding  to  the  low  and  high  ray  modes  of  reflection  from 
the  F-layer.  We  attempted  to  reproduce  these  observed  widths  by  varying 
the  parameters  of  the  DNA  scintillation  model  in  a  fit  of  the  data. 
Initially  we  fit  the  low  ray  modes  exactly  by  varying  the  spectral  index 
and  plasma  velocity  parameters,  fixing  the  other  model  parameters  to 
nominal  values.  When  these  resulting  parameters  were  applied  to  the  high 
ray  modes,  predicted  delay  widths  and  doppler  spreads  were  consistently 
too  large.  It  was  observed  that  in  all  but  one  case  these  spreads  could 
be  brought  back  into  substantial  agreement  with  the  data  by  allowing  one 
parameter,  the  fractional  electron  density  variance  £  =  <  ANp  ,  to 

take  on  different  values  for  the  high  and  low  rays.  This  result  was 
assumed  to  be  an  approximation  to  a  single  spatially  dependent  £.  When 
the  high  and  low  ray  modes  were  simultaneously  fit  by  varying  the  spectral 
index,  plasma  velocity,  and  fractional  electron  densities  for  the  high  and 
low  rays,  the  results  were  found  to  be  in  good  agreement  with  the  data. 
Furthermore,  the  resulting  spectral  indices  were  very  near  those  inferred 
from  the  DNA  HILAT  measurements  carried  out  in  the  same  region. 


We  have  scaled  the  oblique  ionograms  of  the  HF  channel  Probe  by 
fitting  them  directly  to  the  delays  produced  by  a  quasi -par abol ic  layer. 
This  method  avoids  the  approximate  step  of  transforming  the  oblique  iono- 
gram  to  its  vertical  equivalent.  The  three  parameters  of  the  quasi¬ 
parabolic  layer;  the  critical  frequency,  height  of  layer  maximum,  and 
layer  semi-thickness  were  supplemented  with  a  fourth  parameter,  an  abso¬ 
lute  delay  shift  of  the  ionogram  data.  This  was  done  to  account  for  a  200 
usee  uncertainty  in  the  absolute  delay  measurements  of  the  HF  Channel 
probe.  We  found  that  fits  to  the  layer  obtained  by  including  this  fourth 
parameter  resulted  in  better  positioning  of  the  E-layer  trace  when  that 
was  visible  in  the  data.  Still,  these  delay  uncertainties  are  a  trouble¬ 
some  point  in  our  analysis  because  of  the  sensitivity  of  the  DNA  scintil¬ 
lation  model  calculations  to  ray  path  and  plasma  frequency.  For  instance, 
the  success  attained  by  considering  variable  fractional  electron  density 
variances  could  conceivably  be  explained  by  some  consistent  problem  with 
our  model  ionosphere  and  the  subsequent  high  rays  obtained  from  it.  In 
future  campaigns  of  the  HF  Channel  Probe,  this  delay  uncertainty  will  be 
reduced  by  an  order  of  magnitude,  so  this  question  should  be  answered  by 
future  analysis  of  this  new  data. 

Our  investigation  of  the  parameter  dependence  of  the  DNA  scin¬ 
tillation  model  showed  that  of  the  six  model  parameters,  the  spectral 
index  is  substantially  the  most  sensitive.  Order  of  magnitude  changes  of 
the  calculated  spreads  occurred  as  the  spectral  index  varied  over  the 
range  of  values  inferred  from  the  DNA  HILAT  measurements.  The  model  show¬ 
ed  less  sensitivity  to  the  other  parameters.  The  perpendicular  structure 
outer  scale  size  L  ,  for  instance,  produced  variations  in  delay  and  dopp- 
ler  spreads  of  less  than  a  factor  of  two  when  ranqed  over  reasonable 
values.  The  plasma  velocity  parameter  only  enters  the  doppler  spread 
calculation  and  has  a  linear  dependence.  The  delay  spread  was  shown  to 
depend  on  the  fractional  electron  density  variance  £  in  a  nearly  linear 
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fashion.  The  doppler  spread  scaled  essentially  as  the  square  root  of  this 
parameter.  The  remaining  parameters,  the  parallel  structure  outer  scale 
size  L  and  the  ratio  of  inner  to  outer  structure  scale  size  e,  were  shown 
to  have  little  effect  on  the  calculated  spreads  for  any  reasonable  range 
of  values. 

Our  initial  fit  of  the  low  ray  spreads  by  varying  only  the  spec¬ 
tral  index  and  plasma  velocity  and  the  subsequent  application  of  these 
parameters  to  the  high  ray  yielded  results  that,  even  though  the  hiqh  ray 
spreads  were  consistently  too  large,  were  accurate  to  an  order  of  magni¬ 
tude  or  so.  This  provides  some  validity  to  the  scintillation  model  as  a 
predictive  tool  for  highly  oblique  paths  at  HF.  By  fixing  most  of  the 

model  parameters  to  nominal  values  and  varying  the  spectral  index  over  a 

sensible  range,  say  SI  *  1.5  to  2.5,  expected  delay  spread  bounds  can  be 

reasonably  estimated,  especially  for  low  ray  modes  which  are  the  predomi¬ 
nant  modes  of  interest  for  HF  communication. 

Interpreting  the  good  agreement  with  the  data  obtained  in  ou*- 
four  parameter  fit  summarized  in  Table  10  is  more  speculative.  As  we  have 
already  remarked,  agreement  might  also  have  been  obtained  by  alter  do 
layer  parameters  to  vary  the  high  ray  path  or  by  choosinq  r  >  •• 

spectral  index  or  Lx  spatially  dependent.  It  is  important  *  ■  -  •  * 

ever,  that  by  making  this  single  parameter  r,  effectively  • 

dent,  both  the  calculated  high  ray  delay  and  Hopple1-  ■  . 

c  i  led  to  the  data.  Also,  freeing  up  as  a  fit  pa*-  r>-  •  •  • 
fit  spectral  indices  which  are  very  dose  »n  *n.  v  ■ 

As  yet  we  have  no  theoretical  iust  i  f  i  <  at  '  er  *  • 
density  variances  should  he  relatively  ■ 

rouqh  test  where  we  examined  the  e\  ••  • 
low  and  high  rays  for  earh  r  io-,  *.  • 

the  magnitude  of  the  qrad'en*  *  •>.. 


six  data  sets  could  be  explained.  Theoretical  work  and  further  examina¬ 
tion  of  the  data  is  obviously  required  to  substantiate  any  of  these  specu¬ 
lations. 

The  DNA  scintillation  model  has  been  shown  to  be  capable  of 
reproducing  the  spreads  in  delay  and  doppler  observed  by  the  HF  Channel 
Probe  in  the  polar  region  for  quite  reasonable  values  of  the  model  param¬ 
eters.  It  must  be  born  in  mind  that  these  tests  were  for  highly  oblique 
HF  propagation  paths  only  (i.e.,  launch  elevation  angles  less  than  twenty 
degrees  or  so),  and  no  test  of  the  model  has  been  made  in  situations  where 
the  local  plasma  frequency  of  the  ionosphere  is  close  to  the  wave  carrier 
frequency,  as  would  be  the  case  in  nearly  vertical  paths.  The  model,  how¬ 
ever,  showed  no  gross  failure  when  applied  to  oblique  HF  propagation. 
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